PM/00-41 
November 2000 



Three— Body Decays of Top and Bottom Squarks 



A. Djouadi 1 and Y. Mambrini 2 



Laboratoire de Physique Mathematique et Theorique, UMR5825-CNRS, 
Universite de Montpellier II, F-34095 Montpellier Cedex 5, France. 

1 djouadi@lpm.univ-montp2.fr 
2 yann@lpm.univ-montp2.fr 



Abstract 

We investigate the decays of third generation scalar quarks in the Minimal Super- 
symmetric extension of the Standard Model, focusing on the three-body modes. 
We calculate the partial widths of the decays of heavier top and bottom squarks 
into the lighter ones and a fermion pair [through virtual vector boson, Higgs boson 
or gaugino exchanges] and the partial widths of the three-body decays of both top 
squarks into bottom quarks and a pair of fermion and scalar fermion [we consider 
the case of lighter f or b states] and into a bottom quark, the lightest neutralino 
and a W or a charged Higgs boson H^. Some of these decay modes are shown to 
have substantial branching ratios in some areas of the parameter space. 
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1. Introduction 



In the Minimal Supersymmetric extension of the Standard Model (MSSM) [0, the spin- 
zero partners of third generation standard chiral fermions can be significantly lighter 
than the corresponding scalar partners of first and second generation fermions. This is 
essentially due to the relatively large values of third generation fermion Yukawa couplings 
which enter in the non-diagonal entries of the sfermion mass matrices, the diagonalization 
of which turn the left- and right-handed current eigenstates and into the mass 
eigenstates f\ and fi . The mixing can generate a sizeable splitting between the masses 
of the two physical states and leads to a lighter sfermion f\ with a mass possibly much 
smaller that the masses of the other sfermions. The situation can be even more special 
in the case of the lightest top squark, t±, whose mass can be smaller than the one of its 
partner the top quark, m^ < m t , to be compared with the experimental lower bound on 
the masses of the first and second generation squarks, rriq > (9(250 GeV) 0. 

The fact that the top quark is heavy leads to distinct phenomenological features for 
the decays of its scalar partners. Indeed, while the other squarks can decay directly into 
(almost) massless quarks and the lightest neutralino x?, which is always kinematically ac- 
cessible since in the MSSM the neutralino Xi is assumed to be the lightest supersymmetric 
particle (LSP), the decay channels ti —>■ t\\ are kinematically closed for mi. < m t + m x o. 
If, in addition, m~ ti < m& + m x + with xf being the lightest chargino, the decay modes 
U ~ > bxf are not accessible and the only two-body decay channel which would be allowed 
is the loop induced and flavor changing decay into a charm quark and the LSP, ti cx\ 
[[|. The other possible mode is the four-body decay channel into a bottom quark, the 
LSP and two massless fermions, ti — > bxiff, which occur through virtual exchange of 
top quarks, charginos and scalar fermions |5[]. 

For relatively heavier top squarks, the three-body decay channels 

U - bW + x°i , bH + Xi (1) 

where H ± is the MSSM charged Higgs boson, can be accessible; Fig. la-b. These decays 
have been discussed in Ref. || [F] in the case of the lightest top squark and have been 
shown to be [at least for the one with W boson final states] often dominant in the case 
where mi x < mt + m x o and m& + m x + . In addition, if sleptons are lighter than squarks 
[as is often the case in models with a common scalar mass at the GUT scale such as the 
minimal Supergravity (mSUGRA) model] the modesf] 

U — > bl + ui and/or bl + ui (2) 

become possible 0, 0, jl0|l , Fig. lc. In the case of the lightest top squarks, they can be 
largely dominating over the loop induced CX\ mode. 

In this paper, we point out that these three-body decay modes are important not only 
for the lightest top squark, but also for the heavier one. In addition, we investigate a new 

1 This mode has also been discussed in Ref. B for first and second generation slepton decays into 
lighter r sleptons in the context of gauge mediated Supersymmetry breaking models H . 
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possibility which is the decay of the top squarks into a fermion-antifermion pair and the 
lightest b state, which is mediated by the virtual exchange of W and H + bosons: 

U - hff (3) 

bi can become the lightest scalar quark in the case where the ratio of the vacuum expec- 
tation values of the two-Higgs doublet fields in the MSSM, tan j3, is largeQ. 

For the heavier top squark, t 2 , another possibility would be the three-body decay into 
the lightest top squark and a fermion pair [with / ^ b] through the exchange of the Z 
and the MSSM neutral Higgs bosons [the CP-even h, H and the CP-odd A bosons], 

h - tiff (4) 

These modes apply also for the charged decays of heavier bottom squarks into top squarks 
(and vice-versa) which, as in eq. (3), occur through W and H + boson exchanges^ 

h - q'jff (5) 

For bb final states, one needs to include in the case of t% — > t\ bb the contributions of the 
exchange of the two charginos states Xi2i Fig- ld-e. This is also the case of the decay 
mode, 62 —>■ b\ bb, where one has in addition, the virtual exchange of neutralinos and 
gluinos, Fig. le, which have to be taken into account. The latter process is a general- 
ization [since the mixing pattern is more complicated] of the decay modes of first and 
second generation squarks into light scalar bottoms discussed in Ref. ||l4jl , and would be 



in competition with at least the two-body mode b 2 — > The latter channel is always 
open since Xi is the LSP, but the b-b 2 -Xx coupling can be small, leaving the possibility 
to the three-body mode to occur at a sizeable rate. 

In this paper we analyze all the three-body decay modes, eqs. (2-5), discussed above. 
We will give complete analytical expressions for the Dalitz plot densities in terms of the 
energies of the final fermions as well as the fully integrated partial decay widths. In 
addition, we investigate the t\ and t 2 decay modes of eq. (1) which have been already 
discussed in Ref. [0 for the lightest top squark t\. In this case however, only the Dalitz 



densities will be given; the more complete and lengthy formulae can be found in Ref. [15 

The rest of the paper is organized as follows. In the next section, we will discuss the 
main properties of top and bottom squarks and summarize their two-body decay modes 
for completeness. In sections 3 and 4, we analyze respectively, the decay modes of top 
and bottom squarks into lighter sfermions and fermion pairs, and the decays of the two 
top squarks into neutralinos, b quarks and W or H + bosons. A numerical illustration is 
given in section 5 and a brief conclusion in section 6. 



2 The scenario with large values of tan/3, tan/3 ~ rnt/rrib, is favored in models with Yukawa coupling 
unification at the GUT scale jy]; the other possible solution, with tan/3 ~ 1.5, seems to be ruled out by 
the negative searches of MSSM Higgs bosons at LEP2 E^j. 

3 If the mass splitting between the initial and final scalar eigenstates is large enough, the gauge and 
Higgs bosons become real, and we have the two-body decays into gauge and Higgs bosons which have 
been recently analyzed in Ref. tt3]. 



3 




d) 




Figure 1: Feynman diagrams for the three-body decay modes of top and bottom squarks. 
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2. The Two-Body Decay Modes 



In this section, we will summarize for completeness the two-body decays of scalar quarks. 
This will give us the opportunity to exhibit the various couplings of squarks to charginos, 
neutralinos, Higgs and gauge bosons which will be needed later on, and to discuss the 
third generation sfermion mass spectrum and mixing pattern. 

2.1 Sfermion masses and mixing 

As mentioned earlier, the left-handed and right-handed sfermions of the third generation 
fi and fn [the current eigenstates] can strongly mix to form the mass eigenstates f\ and 
$2'-, the mass matrices which determine the mixing are given by 



Ml = 



m LL m fA f 
m,A f m\ R 



(6) 



with, in terms of the soft SUSY-breaking scalar masses mj L and mj R , the trilinear cou- 
pling Af, the higgsino mass parameter /i and tan/3 = vu/vd, the ratio of the vacuum 
expectation values of the two-Higgs doublet fields 



m 2 



ll = m} + mj L + (lj-e f s 2 w )cos2(3M 2 z 

2 , 2 
f +m fn 



m RR = rn 2 f + m 2 f + e s s 2 w cos 2/3 Mf 



A f = A f -^(tan/3)" 27 / (7) 

with ef and If the electric charge and weak isospin of the sfermion / and sly = 1 — = 
sin 2 6w- In the stop sector, the mixing is strong for large values of the trilinear coupling 
A t and/or for large values of /x with small values of tan (3. In the case of the scalar bottom 
and tau lepton, the mixing is large when tan (3 and the parameter \i are large. 

The mass matrices eq. (6) are diagonalized by 2 x 2 rotation matrices of angle Of 





VJ = \ 9f 6f , c 8f = cos 9 f and s 9 =sin9 f (8) 
The mixing angle Of and the sfermion eigenstate masses are then given by 



2.2 Two— body decays into neutralinos and charginos 

If the scalar quarks qi are heavy enough, their main decay modes will be into their partner 
quarks and neutralinos, q~i — > qy^ [j=l-4], and quarks and charginos, q\ — > </Xj [j = 1— 2] . 
The partial decay widths are given at the tree-level byQ: 



q'xt) 



4m? 



aA2(m?,m2,,m2 + ) „ 2 „ 2 

— (4 + 4? )( m « ~~ "V ~ m x +) _ 4 44/ TO 9 ,T7 V 



4mi 



[ID 



where A(x, ?/, 2) = x 2 + y 2 + z 2 — 2 [xy + xz + yz) is the usual two-body phase space 
function and e Xj is the sign of the eigenvalue of the neutralino x?- The couplings a i3 - and 
by for the neutral decay are given by 



°2j 



m / r / J S <>, 1 _ J 



V2M w s w \ c 6f 



in 



fj_\ 



c'e, 



V2M w s w \ sof 
with r t = Nj±/ sin (3 and r T = r b = Nj 3 / cos (3 and 



— e 



CO, 
-8e f 



:i2i 



-Rj 



V2 



N' n + (if - e f s 2 w 



c w s w 



-V2 



N' 3l - ^ N> 2 

J c w J 



(13) 



while the couplings Cj,- and gL- for the charged decay mode are given, for U decays, by: 



and 



^3 

J 
Xi 

it 

l 2j 



m b Up 



C0i 



\/2M w s w cos/? \ SQ t 



v n \ -co* 



mtVj 



J'2 



V2M w s w sin (3 \ ce t 



(14) 



d~ b 
a lj 

d b 

a 2j 



m t Vi 



.ft 



V2M w s w sin 13 \ se b 



Uji 
s w 



~ c e b 
s e b 



+ 



m b U j2 



\/2M w sw cos/3 



s e b 



J The QCD corrections to these decay modes have been calculated in Ref. faj. 



(15) 
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for hi state decays [in the case of r sleptons, one has to replace in the previous equations, b 
by r and to set m t = and 6 t = 0]. In these equations, N and U /V are the diagonalizing 
matrices for the neutralino and chargino states with 



N' n = c w Nji + s w N j2 



Sy/Nji + CwN j2 



(16) 



In the case of top squarks, these decays might be not accessible kinematically, and the 
only allowed two-body decay will be the loop induced and flavor changing decay mode 
into a charm quark and the lightest neutralino, £j — > c%°. To a good approximation, the 
partial decay widths [in mSUGRA] are given by [|J: 



r& - ex") 




J Lj\ 



2M 2 v s 2 w cos 2 (3 



log 



^GUT ] 
v M W J 



A; 



X 



2 2 

m~ — m~ 



(17) 



Ai = -c et {m\ L + m\ R + m 2 Hl + Af) + s dt m t A b 

A 2 = s et (m 2 L + m\ R + m 2 Hi + A 2 b ) + c 6t m t A b (18) 

The widths are suppressed by the CKM matrix element V c b ~ 0.05 and the (running) b 
quark mass squared ml ~ (3 GeV) 2 , but very strongly enhanced by the term log (Aq UT /M^) 
with Aqut — 2 • 10 16 GeV. Assuming proper electroweak symmetry breaking, the Higgs 
scalar mass can be written in terms of //, tan (3 and the pseudoscalar Higgs boson 
mass M A as m 2 Hl = M\ sin 2 (3 - cos 2(3 - fi 2 . 



2.3 Two— body decays of q 2 into gauge and Higgs bosons 



If the mass splitting between two squarks of the same generation is large enough, the 
heavier squark can decay into a lighter one plus a gauge boson V = W, Z or a Higgs 
boson $ = h,H, A, H ± . The partial decay widths are given at the tree-level by0: 



a 



-, # A 1/2 (m|, M|,m|) 



4m|M 2 
a 2 



(19) 
(20) 



In these equations, the couplings of the Higgs bosons to squarks, g^q'.®, read in the case 
of neutral Higgs bosons: 



9qiq~2h 



As w M w . 



M 2 z s 26q (2/J - Ae q s 2 w ) sin(a + (3) + 2m q c 2Qq (A q r q 2 + 2/J /x r 



(21) 



9q\q~2H 



AswM w 



M 2 z s 2Q (2ll - Ae q s 2 w ) cos(« + /3) + 2m q c 29 (A q rf - 2I 3 ur q 2 ) 



(22) 



3 The QCD corrections to these decay modes have also been calculated and can be found in Ref. 1 18 
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9q\qiA 9qiq\A 



2s w M w l 



(23) 



with the coefficients r\ 2 as [a is a mixing angle in the CP-even Higgs sector of the MSSM, 
and at the tree-level, can be expressed only in terms of Ma and tan (3] 



t sin a 



cos a h cos a h sin a 
*i = a , ^2 _ 



sin/5 sin/3 cos/3 cos/3 

In the case of the charged Higgs boson, the couplings to squarks are given by 

( 2 



2s w M w ^ V Jik <* H \ Jij 



(24) 



(25) 



with the matrix C^ H ± summarizing the couplings of the H + bosons to the squark current 
eigenstates; it is given by 

„ n: ( m\ tan/5 + m 2 /t&n[3 - M^ sin2/3 m b (A b tan(3 + //) \ , . 

° to± : V ^ ^ mj ( A/ tan (3 + fi) 2m t m b / sin 2(3 ) [ b) 



Turning to the couplings of squarks to the W and Z gauge bosons, one has 

QqifaZ — -9q 2 qiZ ~ ~ — \*t) 



4s w c w 



1 



9q&W 



c dq c e > 



-ce q s e > q 



V2s w V - s e q °e' q se q s e > q 



(28) 



Finally, for the next section, we will need the couplings of the W, Z gauge bosons and the 
four Higgs bosons h, H, A and to fermions [r{ 2 are defined above]: 



v ffz 



9ffh 



2Ij-4e f s 2 w 

4:C W S W 

f 

rrif r J 2 



a ffz 



2s w M w ' 9ffH 2s w M w ' JJJ " 2s w M w 
s ma tan f3 + m u cot(3 p m^tan/S — m u cot/3 



2/) 



4c w s w 
f 

mfr{ 



V ffW - a,ff W 



2y/2s w 



9ffA 



m/(tan/3) 21 f 



9udH± — 



2V2s w M w ' 9udH± 2V2s w M w 
and the couplings of W and H + bosons to chargino/neutralino pairs: 



9 v o 



-^-[-N^ + V^N^] ,g* 



(29) 

(30) 
(31) 



xgxjw- 



-J— [N: 3 U J2 + V2N* 2 U n } (32) 
y2s w 



9 v o 



n R 



x°xJ H ~ 



cos (3 

sw 
sin (3 

sw 



NisUj! - -j= (N i2 + tan^^i) U j2 



(33) 
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3. Decays into scalar fermion final states 

In this section, we will analyze the decay modes of top and bottom squarks into lighter 
scalar fermions, their quark partners and bottom quarks. We will neglect for simplicity 
the masses of the final state fermions [except in the Yukawa couplings] since, even in the 
case of the bottom quark and tau lepton, it is a very good approximation for the t and 
b masses of the order (9(100 GeV) that we are considering. In the case of top quark 
final states, the t-mass effects have of course to be taken into account and they can be 
found in Ref. fI5 |. However, in this case and for top squarks for instance, the two-body 



decays ti — > txi are kinematically allowed and will dominate over all other decays. Thus, 
throughout this paper, we will not consider t-quark final states and take all the fermions 
to be massless. 

3.1 Scalar top decays into lighter sleptons 

We start by considering the decay of top squarks^ into lighter sleptons, U — > bllj with lj 
being either a sneutrino lj = i>i or a charged slepton eigenstate ii )2 [I = e, fi, r], which 
occurs through the exchange of the two chargino states, Xi2- m terms of the reduced 
energies of the final state particles and the reduced slepton mass defined as: 

xx = H p l% f b , x 2 = 2 Pii Pl , x 3 = 2 Pil = 2 — xi — x 2 , /ir = — % (34) 
m~ m~ m~ m~ 

Li Li l<i l<i 

the Dalitz density of the decay mode reads: 

^Bl^ mh t ,35, 

dXidX2 l07T £jz x 

z lk C il cl jk cl jl + ^ik^il^jk^jl) dG[ kl + y/ ^Xk^Xi ( c ifc c ^jffc^j7 + ^ik^il^jk^jl) ^2kl 

with the two functions dG\ kl and dG l 2 ki are given byQ 

dG f = (l-x 1 )(l-x 2 )- 

dGtij — -jz. rjz ^ r (36) 

(1 -Xx -fi Xi )[l -x x -fin) 



6 The expressions that we will write in this section are also valid, with the proper change of the 
couplings and masses, in the case of bottom squark decays which occur through neutralino exchange. 
This decay is of importance in models where SUSY is broken by gauge interactions (the so-called GMSB 
models Q), and where the lightest SUSY particle, the gravitino, couples very weakly to matter. The b 
states will then mainly decay through this channel (but with neutralino exchange) into the tau slepton 
which is in general the next -to-lightest SUSY particle; see Ref. Q. 

7 These functions are the same as the ones appearing in the simpler case of first and second generation 
squark decays given in eq. (5) of Ref. Q. Note that there is a typographical error in the first part of 
the latter equation: fj,^ has to be replaced by — /ig; the integrated form, eq. (8) of Ref. G4j, based on the 
correct Dalitz density is the same as eq. (38) of the present paper. 
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Integrating the functions over the phase space, with boundary conditions: 



1 — Xi — fif < x 2 < 1 



1 — X\ 



, < Xi < 1 — flf 



one obtains the functions G U j and G 2i j which read: 



(37) 



rf 



- (/i / -l)(3/i / + 3-2/i Xi -2/i Xj )-2 



A%A% 



■log (if 



(A*/ - Vxi) 2 (v Xi ~ l ? Vxi - A*/ _ 2 W - VxjfiVxi ~ l ? A^x, - A*/ 



(38) 



log 



A% - 1 



log 



A*x, - 1 



Vxi^Xj 

vxMf-vxiYiv 



A*/(A*xi + A* 



2 A%A^i 



log/// + -(1 - Hj) (jif + VxiVxj) 



(39) 



l ) 2 , A% - A*/ A%(a*/~ - IbaYUha ~ l f , A^ ~ A*/' 
— log ^ p — ^— log ■ 



2 A%(a% - A%) A% " 1 2 A%(a^ - A%) ~~° Iba " 1 

In the case where Xj — Xi — X [i- e - for the squared terms], the expressions simplify to: 



GL = ^ / -l)(5-6^ + 5|i / -2^)-^log|i / 



+ 2 



4^/ 
(/y ~ A^x) 



(A*x - 1 )(A*/ + A*/A*x + A^x - 3 a4) log 



A^x 



x 



(40) 



r f _ 



+ 



Tl ~ ^x - 2 a4 - 2 A*/ + A*/A*x) + A*/(l + A*/ - — ) logA*/ 



ACI 2 



A^x 



(a*/-A*x)(a*x- ^(A*/-^) , A* x — 1 



A^x 



log- 



(41) 



3.2 Scalar top and bottom decays into lighter squarks and fermion pairs f ^ b 

The neutral decays q 2 — > with q = t or b, are mediated only by Z boson and Ai, iJ, A 

boson exchanges if the final state fermion is not a partner of the decaying squark. The 
Dalitz density, with the reduced energies x\ = 2(pq 2 ■Pf)/ m 1 2 and x 2 = 2(Pq 2 'Pf)l 
and the reduced mass squared jiq = rn^Jm 2 ^ is given by: 

dr(g 2 - qjf) a 2 N c 



m 



92 



dxidx 2 



8vr 



51 9qiq2®9qiq2®'9ff®9ff®' d-F|$/ 
L &,Q'=h,H 

+9lq 2 A9 2 ff AdFf A + ^lq 2Z (v 2 ff z + a)jz) H? z 



(42) 



For the charged decay mode, qi — * q'jff [i-e. those of the decays ti <2 — > b\ 2 ff' and 
6i ;2 — > tipff which are allowed by phase space], mediated by W and H + boson [if the 
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final fermion / is not a partner of the squark q'j], the Dalitz density reads 



rrin 



+9 



4#- ~, nrWfW + « 



$ i q' i W\ u ffW ^ u SfW 
2 i I „P \2 



wiy 



dF H±H± 



(43) 



where xi and X2 are as above and the reduced mass is now fiq = m 2 ~, /ni-.; Nc is the color 
factor of the fermion /, Nq = 3 for quarks and Nq = 1 for leptons. 

The two functions for the exchange of gauge bosons and scalar bosons dFy V [V = 
Z, W] and dF| $ [$, $' = h, H, A, H±] are given by: 

(1 - Xl)(l - X 2 ) - Hq 



dFL, 



(Xt + X 2 - 1 + /ig — Hv) 2 

X\ + X 2 - 1 + jlq 



(44) 



(45) 



{X\ + X 2 — 1 + fig — /i<J.)(Xl + X 2 — 1 + /ig — /!$') 

Integrating over the phase with boundary conditions as in eq. (37), and using the phase 
space function, 

\{lix, Hy) = -1 + 2/ix + 2/iy - (fix - /iy) 2 (46) 

one obtains the integrated functions^ for the partial decay widths [which have to be 
multiplied by the same factors as in eqs. (42,43)] 



5(1 + llq) ~ 4/iy H \(fl V , fiq) +[\(PV, Pq) ~ 2 Vq ) fog Pq 

//y J 



+2(1 - nv + /Xg)\/A(/xy,/Xg)Arctan 



/iy(l - /iy + /ig) - A(/iy, fig) 



(47) 



(/ig - 1) + -(1 + fig ~ fi* ~ HQ,) \Og llq 



Arctan 



1 - Hq)\JK^,Hq) 



~\~ Arctan 



/i<j>(l - yU* + A«g) - A(/X$, /ig) 

(1 



(48) 



The latter function reduces in the case where $ = i.e. for the squared terms, to 



2(/iq -!) + -(! + /ig- 2/1$) log ^ 



/i$(l - fl<j> + Hq) + \{flq,fl<S>) 

A(/i«j, 



Arctan 



1 - ^V^^'ft) 



./i$(l - )U$ + /iq) - A(/Z$, //g) 



(49) 



8 Note that the function -Fyy is the same as the one obtained in Ref. [jl9| for the three-body decays 
of a heavy Higgs boson into a lighter Higgs boson and a fermion-antifcrmion pair. 
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3.3 Scalar top and bottom decays into lighter squarks and bb pairs 

In the case of the decays t 2 ~~ * t\bb, there are additional contributions with the exchange 
of charginos, while in the case of the decay b 2 — * b\bb one has to include the contribu- 
tions of virtual neutralinos and gluinos. The Dalitz density eqs. (42,43) have then to be 
transformed according to: 



dr(t~ 2 -> hbb) dT(q 2 



da^da^ 



dx x dx 2 



qiff) f=b ^ o? 

l07T 



q=t 



2 

E 

k,l=l 



i c \k c \l c 2k c 2l + f 4fc c 4z f 4fc C 4j 



ti 
lkl 



i C lk C lld t 2kd t 2l + ^lfc^U C 2fc C 2z) dGJkl 



-4E 



k=l L 



9tihZ [ C lk C lk( V bbZ - Ctbbz) + $2k$\k{ v bbZ + a>bbz)) dGy fc 



-2 ^2 dbb^giih^lk^lk + d 



2k C lk, 



H x + dG 



(50) 



dr(6 2 - b x bb) dY{q 2 - ft//) /=b of 

3tt m&2 



dxidx 2 



dxidx 2 



g=6 



2s i , c ! , ( lGi^5 + (ci + si )^gdG25i 



16tt 



6 2 



E 

kJ, 



a \k a \l a 2k a 2l + ^lfc^lZ^2A;^2/) dG^jy + {P'\k a \^ > 2k\l + ^\k^\l a 2k a 2l) 



/x x o/x x o dG 2 1 fc/ 



-4E 



c 2k C lk( V bbZ — Q-bbz) + d 2 kd\ k (VbbZ + d-bbZ 



b jb 



dG bl 



Vk 



"2E5 ,b6 *^6i62*( C 2fc^lfc + d 2 k c lk) \/Jh^ k dG^ 



(51) 



Note that the sums run only on the virtual states; for instance, in the case of the decay 
62 — » one has to discard the exchange of the lightest neutralino Xi [since it is the 
LSP and the decay 62 — > bxl always occurs at the two-body level if ra& = 0] and add the 
two-body partial width r(6i — > bxl) to the total decay widthP]. Note also that the gluino 
exchange diagram does not interfere with the other diagrams due to color conservation. 

The functions dG{ kl and dG^i have been given previously, while the new functions 



dG^j and dG^ are given by: 



dG q Vi 



(1 - Xi)(l - X 2 ) - Hq 



dG q 



{x l + X 2 - 1 + fiq - yUy)(l ~ X X - fl Xi ) 

X 1 +X 2 — l + fiq 
(X 1 +X 2 — 1 + Hq ~ yU$)(l -Xi- fJL Xi ) 



(52) 



9 If the finite widths of the exchanged particles are consistently included in the expressions, one can 
use them also for on-shell exchanged particles. However, in the case of the decay b\ — > bxi this procedure 



has always to be done since \i is stable. 
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When integrating over the phase space, one obtains: 



0« 



-(1 + piq + 2/iy) + (1 + /J,q - 2/i Xi - 2/iy) log /iy 



+ log 



- Hxi)(Vxi - 1 ) ~ ^[K^v, Hq) + 2 ^ Xl (l - //v + /ig) - 6/i g -] log/^~ 



+ nW'W 1 + ^" 2 ' 1 * -//y Arctan — — ? r 

2 V \^A(/iv,//g) +^y(/iy - 1 - // $ J 

V „V\ 



(53) 



G li = (N - 1) - — log/Xqf + 



log- 



+ M*/(^,r*r?) (54) 



with the function /, with arguments r± = | 1 + /^x — A 4 ? + J—\(/ix, A*g) > defined as: 



/>, it, u) = -/(l) + f(u) + /(«) + log 3 log 



(55) 



A i x» - 1 

where in terms of the Spence functions defined as, Li 2 (x) = J x dtt -1 log(l — t), one has: 



/(*) = log (^ - 1 + x) log -i-^i- + Li 2 ( ^-g 



Lio 



(56) 



4. Scalar top decays into W, /f + bosons and the LSP 

In this section we will analyze the three-body decay modes of top squarks, tj — > &X?-B 
with B = W, H + . This is a (straightforward) generalization of the modes t\ — > 
discussed in Refs. 0, [5J since here, we will consider both top squarks in the initial state 
and any neutralino in the final state. Here again, we will neglect the 6-quark mass in the 
amplitude squared and in the phase space, as well as the finite widths of the exchanged 
particles [the latter can be easily included in the propagators]; the complete expressions 
with a finite m& value in the propagators and in the phase space [which gives a better 
approximation for stop masses of order 100 GeV] can be found in Ref. [|l5|] . 

In terms of the reduced energies of the final particles X\ = 2(p: ■ pb)/mf , x 2 = 
(Pit ' Px])/ m i an d x 3 — (Pii ' Pb)/™^., and for the reduced masses fix = Px/ m % [ we w iU 
drop the index for Xp A*x = /V']' anc ^ introducing the new scaled variables: 



Pb ■ P Xj 



ra- 



in 



Pb-PB 



■y.3 



Pxj ■ Pb 



mt 



(57) 



the Dalitz densities for the decay modes ti —>■ bx®B are given by: 

= ^- [ r g + r* +r* + 2if + 2it + 2 R 

dxidx 2 16tt L bb xx « b x u 

where the terms correspond to the square of the contributions of the sbottom, chargino 
and top quark exchange diagrams and the interference terms. 



(5f 
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4.1 The decay i t -> bx°jW 

In the case of the decay U — > bx°jW, one has for the various terms: 



6b 



k,l=l 



kl 



(59) 



xx 



= E 



a,z=i (l-xi-^+)(l-a;i-^ x +) 1 

42/3(2/1 + 2/2 + vwvm) + y\{vx - m) + 2^*2/2(1 - ^2/3) 



/nrl+ 



+ 



(60) 



3+ 



4j/2j/i(A*wJ/2 + 1) - + 4 2/22/3 



(1 -x 2 +/i x -/i t ) 2 r* j 

+ + %22/3/V 1 ] - 4 v //i t /x x Cjy 2 [3 + 2^2/2] 



(61) 



-nVK 

i 5x 



-4 



y 



X; 



(2/2 + y^ivxVi ~ 22/12/3 W 



+ 2/1 (22/i + 2/2 - 2/3 + /i x ) + VxVi 



+ C k^xJxJ ^ ~ ^wVi{V2 + 2/3) 



(62) 



1 6t 



E 



; (1 - x 3 + n x - m Vi -x 2 + n x - tn) 



2/1-2/2^^(1/2 + 2/3) 



+ Cl + 2/12/2 ( ~ 1 + 2(2/2 + 2/3 W) + Mx2/2 - 2/12/3 - 2y 2 | 



(63) 



= E 



^ (1 -X! -/i x +)(l -X 2 



2/22/3) 



+ eg 



8+ 



2/1(22/3 + 22/2 + 42/i - ^w) + 2/2(42/3 + M x ) - 2^2/2(22/12/3 - ^2/2) 



- Zy/mHxC^{y\ + 2/2) - ^jUxUx+Cl 2/2(3 + 2^2/2) 
The various combinations of couplings C ..g read as follows: 
C/°fc = 9t i b k w9t i b l w( ab kj ab ij + bkjtfj) 



(64) 
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r 2 

3± 
Ik 



Jt At /~iL r^L , t t /~fR (~iR 
u ik u il (j jkW (j jlW c ik c il Ur jkW^ r jlW 

it Jt r^R riR , t t (~iL /~iL 
a ik a il (j jkW^ jlW c ik c il (j jkW^'jlW 



rt it (~<L (~<t 



R 



_i_ r t „t nt 
Hk^l^jlW^jkW "T" ^k^l^jlW^jkW 



( a ij) 2 ( a ffw ± Vf f wf + (bifiaffw =f v ffW f 



U lk 

cl 



cl + 



cl- 
ef 
cl- 
ef 
ct 
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a h°lA a2 ffw + v 2 ffW ) 



^kj^ilCjiw + ^kj C ilCflw)9 



(tftjCaCfiw + ab kjdiiGf lw )g 



<Ub k W 



jlWMtib k W 



a ij( v ffw + affw) + hj b ij( v ffw - affw) 



9ti'b k W 



9ti~b k w 



a 



ffWj 



'!■/'!,( 'V/ir - a ffw ) + a b k j b U v ffw + a ffw ^ 

a \jd\kG%w( v }}w + a ffw) + b \jc\ k G jkW {vff W - 
at ij c lkCf k w( v ffw - a ffw) + b%d\ k Gf kw {v f f W + a ffW 
at ij c \kCfkw( v f fw - affw) + %d\kC% w { v f W + affw 
a ijdikCfkw( v ffw + a ffw) + blj4 k Gf h 1 " - " 



'IjdlkCfkwivffw + affw) 
ijCikGfkwivffw - a ffw) 



(65) 



Note that there are two typographical errors^ in the corresponding expressions for 
these amplitudes in terms of the four-momenta, given in Ref. M. In eq. (18) for the 
amplitude square of the chargino exchange contribution, the term 2m 2 should be absent 
in the first square bracket. Furthermore, in the first square bracket of eq. (19) for the 
amplitude squared of the top quark exchange diagram, A/M w x (p x o -pw)(Pb -Pw) 2 should 
be replaced by A/M w x (p x o ■ p b )(p b ■ p w ) 2 , +3M w (p x o ■ p b ) by -M w (p x o ■ p b ) and -A{jp b ■ 
Pw)(p x » ■ Pw) should be replaced by +A(p b ■ pw)(p x <( ' Pw)- [For the numerical analysis, 
the agreement with the figures given in Ref. || is rather good.] 

4.2 The decay U bx]H + 

In the case of the decay ti — > bx°~H + , one has for the various terms: 



1 bb ~ z 2^ U H (-, 

k,l=l v 1 



Hi 



X 3 + fl H - /ig. )(1 ~X 3 + fX H ~ m ) 



(66) 



XX 



fej=i 



(1 -Xi - yU x +)(l - X\ 




2y 3 y2 + yiivx ~ vh) + 2ju x 2/2 



(67) 



v x (yi + yi) [^J\+ D ki + ^ff\+ D 
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We thank Werner Porod for his cooperation in resolving this issue. 
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(1 - x 2 + (I* - fit)' 



D$-(-»hVi + 2y 2 y 3 ) + D\- mi - 2 v ^L>Jy 2 | (68) 



2 D t^K{yi + 2/2) + Dijjr~^ yi 



(69) 



6f 



2E 



^ (1 - x 3 + /x x - /ib fc )(l - x 2 + n x - ih) 



(70) 



1 X* 



E 



^ (1 - x x - /i x +)(l - x 2 + Hx ~ V*) 
+ D s k + (-fi Hyi + 2y 2 y 3 + fi x y 2 ) - ^jMj^D 9 k + ( yi + y 2 ) + ^ X fi x +D 9 k ~y 2 ^j (71) 



with the various combinations of couplings D g given by: 



u m — 

Di + = 

D} k ~ = 

Dl = 

Df k + = 

n 3 - _ 

u ik — 

D ik = 

Dl = 

Dl = 

Dl + = 

Dl~ = 

Dt + = 

Dl~ = 

Dl + = 

Dt = 



c ik c il^ jkH^ jlH 



9t i b k H9t i b l H( a kj a ' b lj + bkjtfj) 

dlkduGjkH^jiH + 
d\kd\iGf kH Gf lH + clkCtfGjkuGjUj 

jt jt (~iL /~iR , t i r-iR (~<L 
Uj ik Uj il Kj jlH u jkH "I" c ik c il u jlH u jkH 

(4) 2 (&+ + 9t b H+? + (bi,) 2 (g? bH+ - 9tbH + ) 



13 

j \2 



(alj) 2 (9t b H + ' 9t bH+ Y + (bl„) 2 (g? bH+ + g? bH+ ) 



(9tbH+ + 9thH+) 2 + (9tbH+ ~ 9tbH+) 2 
{a\jd\iG^ lH + bljC^Gf^g^jj 
(b b kj c iiGjiH + a kjdiiGfi H )g^i kH 

b kj a \j{gtbH+ ~ ga>H+) + b b kjbij(gt b H+ + gtbH+ 

\j a \j(gtbH+ + gtbH+) + ab kjb t ij(g'tbH+ ~ 9tbH+ 

at ij^ikGfkH{9tbH+ + gtbH+) + bijdl k Gj kH (g tbH 
at ijd t ikGfkH(g'tbH+ ~ 9tbH+) + b t ijC t ik Gj kH (gf bH+ + g[ bH +) 
a ijdlkGfkH(gt b H+ - 9tbH+) + b t ij c t ik Gf kH (gf bH+ + 5^+) 
a ij c ikGfkH(gtbH+ + 9tbH+) + b t ij d t ik Gf kH (gf bH + - g[ bH +) 



9tib k H 
9i t b k H 

- ga,H+ ) 



(72) 
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5. Numerical illustrations 



A Fortran code called SDECAY p0| has been developed for the numerical analysis; all the 
partial decay widths and the branching ratios for the two-body and three-body decay 
modes of scalar quarks [as well as the decays of charginos, neutralinos, gluinos, sleptons 
and the four-body decays of the lightest to squark] have been implemented. It has been 
interfaced with the programs SUSPECT |2T[| for the calculation of the supersymmetric 
particle spectrum [including the renormalization group equations for the evolution of the 
SUSY parameters and the implementation of radiative electroweak symmetry breaking] 
and the program HDECAY |22|] for the Higgs boson spectrum and couplings |T7| , |23fl where 
the renormalization improved two-loop radiative corrections in the MSSM Higgs sector 
f24|] and the QCD corrections to the Higgs couplings |2l| have been incorporated. 

We begin our numerical illustration with the decays of the lightest top squark t\. We 
will concentrate on the "unconstrained" MSSM, where for simplicity, we use a common 
soft-SUSY breaking scalar mass rriq for the three generations of squarks and m z ~ for the 
three generations of sleptons, i.e. m~ tL = mi R = m~ bL = m bR = rriq and = nif R = 
vriy L = mj. [We will also assume that the mixing between different generations is absent 
at the tree-level, otherwise the decay mode U — > c\°j would already occur at this stage.] 
The mass splitting between the two mass eigenstates will be then only due to the different 
D-terms of mj L and mj R and to the off-diagonal entries of the sfermion mass matrices. 

The mixing is made strong in the stop sector by taking large values of A t ~ 0(1 TeV); 
in this case the mixing angle is either close to 9 t = n/2 (no mixing) or to ±7r/4 (maximal 
mixing) for respectively small and large values of the entry m t A t compared to the diagonal 
entries of the mass matrix. The mixing is strong in the b and f sectors for large values 
of tan j3 and the parameter /i, almost independently of and A T which will be fixed to 
100 GeV. In the gaugino sector, we will make the usual assumption of the unification of 
the gaugino masses at the GUT scale, leading to the relation Mi = Itan 2 ^^ ~ 

In Fig. 2, we show the branching ratios of the decays of the lightest top squark t\ 
as a function of tan/3 for large values of \i = —750 GeV. This implies that the lightest 
chargino and neutralinos are gaugino like for M-i < 300 GeV, with masses m v + ~ m v o ~ 
2m x o ~ M.2 [with a very small variation with tan /3] . We choose a common squark mass 
rriq of (9(500 GeV) which, for the chosen ii and A t values, leads to a ti with a mass 
between 170 and 250 GeV [depending on the value of tan/3, m$ being smaller for low 
tan/3 values]. In this case, mi x is smaller than mt + m x o and m(, + m x + but possibly larger 
than m x o + M w , m x o + M H + , m~ bi or m fl , allowing to some three-body decay channels to be 
open kinematically. Since these three-body decay modes are of 0(a 2 ), they can compete 
with the t\ — > cxi mode which is of 0(a 3 ) modulo the large logarithm log(AQ UT /M^/). 

In Fig. 2a, the common slepton mass is chosen to be relatively small, mj = 280 GeV, to 
allow for t\ decays into staus and the pseudoscalar A boson mass is taken to be relatively 
large, Ma = 250 GeV, implying a too heavy charged Higgs boson, M H + = J M\ + ~ 
260 GeV, for the decay t\ — > bx\H + to occur, rriq and A t are fixed to 450 GeV and 
1 TeV, respectively, while M2 = 250 GeV. This leads to a scalar fermion spectrum, for 
tan/3 = 5 (45), of m k ~ 170 (230) GeV, m 5l ~ 430 (190) GeV and m fl ~ 270 (140) GeV. 
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Figure 2: The branching ratios for the two-body and three-body decay modes of the 
lightest top squark t\ as a function of tan j3 for /i — —750 GeV. 
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For small values of tan/3, tan/3 < 7, the mixing is too strong in the stop sector and 
ti, being too light to have three-body decays, will mainly decay into c\\ final states. 
For larger values of tan /3, the t\ mass becomes larger than m;, + My/ + m x o an d the 
channel t\ — > bx\W becomes kinematically accessible; it will be largely dominant, with 
a branching ratio above ~ 80%, up to tan/? ~ 30. For tan/? close to the latter value, fi 
becomes relatively light and the decay ti — > frf^ z/ opens up and becomes competitive, the 
branching ratio reaching a maximum at tan/3 ~ 40. For even larger values of tan/?, bi 
becomes also light and the three-body decay t\ — > biff will be the leading decay channel. 
For tan/3 > 50, m^ is larger than Mw + and the two body-decay t\ — > b{W opens 
up and will have a branching ratio close to unity [however, at this stage b\ will eventually 
become lighter than the LSP neutralino]. 

In the scenario of Fig. 2b, the common slepton mass is taken to be larger than previ- 
ously, mi = 500 GeV, leading to heavier sleptons [in particular f's] while the pseudoscalar 
Higgs boson mass is chosen to be smaller, Ma = 100 GeV, leading to a lighter charged 
Higgs boson, M H + ~ 126 GeV, which can thus appear in the decay modes of the t\ state. 
The other parameters are taken to be rriq = 450 GeV and A t = 800 GeV, while M2 
is fixed to 300 GeV. This gives a scalar fermion spectrum, again for tan/? = 5 (45), of 
m- tl ~ 200 (300) GeV, m~ bi ~ 440 (190) GeV and m fl ~ 500 (440) GeV. 

In this scenario, since the t\ mass is slightly larger than previously, the channel t\ — ► 
bx®W is already kinematically open for small values of tan /?, and will be the dominating 
decay mode until the channel t\ — > b\\H + becomes kinematically accessible. The latter 
will be largely dominating for 20 < tan/3 < 40, reaching a branching ratio of ~ 90% for 
tan/? ~ 35, until the opening of the decay channel t\ — > biff which becomes accessible 
for tan/3 ~ 35. This channel becomes then quickly the dominant decay mode of £1. Note 
that for tan (3 > 45, becomes larger than + M w and we have the two-body decay 
mode t\ ^b\W which has a branching ratio very close to unity. 

Let us now turn our attention to the "Constrained MSSM" or minimal Supergravity 
model (mSUGRA) [p(| where the soft SUSY breaking scalar masses, gaugino masses and 
trilinear couplings are universal at the GUT scale; the left- and right-handed sfermion 
masses are then given in terms of the gaugino mass parameter M1/2, the universal scalar 
mass mo, the universal trilinear coupling Aq and tan/3. The soft SUSY breaking scalar 
masses and the trilinear couplings at the low energy scale are given by their Renormal- 
ization Group Equations, the one-loop approximations of which are given for instance 
in Ref. |23], The parameter /x [up to its sign] is fixed by the requirement of proper 



electroweak symmetry breaking. In mSUGRA and in the relatively small tan/? regime, 
due to the running of the (large) top Yukawa coupling, the two top squarks can be much 
lighter than the other squarks, and in contrast with the first two generations one has 
generically a sizeable splitting between mf and m t - at the electroweak scale. Thus, even 
without large mixing, t\ can be much lighter than the other squarks in this scenario. 

In Table 1, we show some of the branching ratios of the lightest top squark for the 
fixed values of the gaugino mass M1/2 = 250 GeV and sign(/i) = — and for several values 
of the scalar mass mo = 100, 150, 200 and 300 and several values of tan /? = 4, 10, 20 and 
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30 [this leads to m x + ~ 2m x o ~ 200 GeV, with a slight dependence on tan/5]. The i\ 
mass is fixed to approximately ~ 200 GeV by varying the trilinear coupling A . One 
sees that the branching ratios for some of the three-body decays [the channels t\ — > bx\W 
and bf\v T \ are sizeable. 

For small values of tan/5 and the chosen values of m , f\ is rather heavy, and the 
only three-body decay channel which is available is t\ — > bWx\ and for m = 150 GeV 
the branching ratio is very close to unity. For larger values of tan/5, fi becomes lighter 
and the phase space for the decay into bWx\ * s suppressed so that only the decay mode 
ti — > bfiu T is largely dominating. For tan/5 > 30, fi becomes too light [with a mass below 
the experimental bound] and the electroweak symmetry breaking does not take place for 
values of m , A t leading to a relatively light stop. A few remarks can be made here: 

- For larger values of My 2 , the top squarks [and all other squarks] become rather 
heavy and the two-body decays into bxt and even to txl are kinematically allowed and 
dominate. For smaller values of Mi/ 2 , the chargino xf becomes too light and again, the 
two-body decay channel ti — > bxt opens up. 

- In the studied examples, the parameter fi is always rather large, > 500 GeV, so 
that all Higgs particles [except for the h boson] are relatively heavy with a mass of 0(\fi\). 
In particular, H + is too heavy for the three-body decay t\ — > bH + Xi to occur. 

- b\ is also rather heavy in the studied scenario. It is only for very large values of 
tan/5 that b\ becomes light enough for the decay t\ — > biff to occur. But in this case, fi 
is even lighter and its mass is smaller than the experimental bound of (9(100 GeV). 



tan /5 


m 


BR(6W X ?) 


BR(WZ) 


BR(c X ?) 


4 


150 


0.993 


4- 10~ 2 


3 • 10- 2 


10 


100 


3- 10- 4 


0.915 


0.085 


20 


250 


0.02 


0.81 


0.17 


30 


300 


0.015 


0.63 


0.355 



Table 1: Some examples of branching ratios for the three-body decays of t\ in the 
mSUGRA model for M 1/2 = 250 GeV, sign(/i) = - and m h ~ 200 GeV. 

We turn now to the decays of the heavier top squark, t 2 - In principle, t 2 should 
have the same decay modes as the lighter ti if the two squarks have approximately the 
same mass [which means that the mixing is not too strong if the left- and right-handed 
soft-SUSY breaking scalar masses, m~ tL and mi L are approximatively the same]. The 
branching ratios would be, however, different because of the different couplings. However, 
if the mass splitting between the two stop eigenstates is sizeable, the additional mode 
h — ► hff through Z and neutral Higgs boson exchanges has to be taken into account. 
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This is shown in Fig. 3 where the decays of t 2 are displayed as a function of tan /3, 
for fi = -350 GeV and M 2 = 310 GeV. We have taken A t = A b = A T = -100 GeV and 
a common slepton mass mj = 200 GeV; in the squark sector, we have used a common 
mass rriq = 400 GeV for the first and second generation squarks but non-universal masses 
mi L = m~ hL = 200 GeV and m^ R = 120 GeV to allow for lighter top squarks with masses 
m- tl ~ 200 and m~ t2 ~ 250 GeV for tan f3 ~ 10 [in this case, m~ t2 is lighter than the chargino 
Xi so that the two-body decay t 2 — > bxt is shut] . 

In this scenario, b\ and the sleptons f\ and v are lighter than t 2 so that the three- 
body decays t 2 biff and t 2 — > bf\v T , bv\l [in the figure we sum the branching ratios for 
all sleptons] are kinematically open; the former decay channel is dominant up to values 
tan/3 ~ 20 where the two-body decay t 2 — > b^W opens up and reaches a branching 
fraction close to unity. The decays into bWx\ and bH + Xi, as well as the loop induced 
decay t 2 — > cxi, are suppressed below the percent level [we cut the branching ratio for 
the decay t 2 — > bWx±, which is mediated by sbottom exchange, when b± becomes on-shell 
since then, we have the decay chain t 2 — > b{W — > biff]. The mass splitting between t 2 
and ti is smaller than the Z and Higgs boson masses and the three-body decay t 2 — > 
occurs at a sizeable rate for small and intermediate values of tan /3, reaching a branching 
ratio of the order of 50% for tan/3 ~ 15. 
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Figure 3: Example of branching ratios for the two-body and three-body decay modes 
of the heaviest top squark t 2 as a function of tan/5 for fi = —350 GeV, M 2 = 310 GeV, 
mj = 200 GeV and trilinear couplings A t = A b = A T = —100 GeV. 
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Finally, for completeness, we have also studied the three-body decays of bottoms 
squarks. In this case, since Xi is the LSP and m b is small, the two-body decay channels 
hi — > bxi are always kinematically open so that three-body decays can be hardly com- 
petitive. The only situation where the latter can be sizeable is when the hbxi couplings 
are very tiny This occurs when the lightest neutralinos are higgsino-like [one has then to 
consider both Xi an d X2 states since in this case, m x o ~ m x o ~ \fi\] so that the coupling 
is suppressed by m^/M^. One also needs rather small tan/3 values not to enhance the 
couplings which are proportional to 1/ cos (3. However, even in this case, the three-body 
decay bi — > t-\W* — > tiff for instance has a small branching ratio: it is only in a rather 
limited range of the MSSM parameter space that it reaches the level of a few percent. 
This is shown in Fig. 4, where for the chosen set of soft-SUSY breaking parameters, 
BR(6i — > tiff) exceeds the percent level only for relatively small values of tan/3. 
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Figure 4: Example of branching ratios for the two-body and three-body decay modes 
of the lightest bottom squark bi as a function of tan/3. The soft-SUSY breaking scalar 
masses are taken to be, mi h = = 190 GeV, m f - fl = 350 GeV and = 300 GeV. 



The smallness of the three-body decay rates is mainly due to the fact that in the 
MSSM one has to take into account the experimental constraints on the squark masses 
and on tan j3 from the negative Higgs boson searches at LEP. In an unconstrained MSSM, 
for instance without the unification of the gaugino masses at the GUT scale [see Ref. [^7 



for examples of models], several constraints on the SUSY parameters can be relaxed 
and some three-body decays might become important. In particular, as it has been 



recently discussed in Ref. ||14|1 , the decay 62 - * bibb through virtual gluino exchange can 
be competitive [since it is a strong interaction process] with the decay channel 62 —>■ bxi 
if the gluino mass is not too large compared to the lightest neutralino masses, as it might 
be the case in the models discussed in Ref. 
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6. Conclusions 

We have performed a comprehensive analysis of the decays of third generation squarks, 
focusing on the three body decay modes. Because of the large value of the top quark 
mass and the possible large mixing in the stop and sbottom sectors which leads to sizeable 
splitting between the masses of the two physical states, the decay pattern of scalar top [and 
to a lesser extent bottom] quarks can be dramatically different from the decay pattern of 
first and second generation squarks, which simply decay into their almost massless partner 
quarks and neutralino or chargino states. 

Several new decay channels, including the cascade decays of heavier squarks into lighter 
ones and fermion-antifermion pairs as well the decays of both top squarks into W, H + 
bosons and the lightest neutralinos or into leptons and lighter sleptons, are possible. In 
some areas of the MSSM parameter space, these additional decay modes can have sizeable 
branching fractions, and they can even be the dominating decay modes. These channels 
need therefore to be taken into account in the search of scalar top and bottom quarks at 
present and future colliders. 
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